Tunneling between Dilute GaAs Hole Layers 
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We report interlayer tunneling measurements between very dilute two-dimensional GaAs hole 
layers. Surprisingly, the shape and temperature-dependence of the tunneling spectrum can be 
explained with a Fermi liquid-based tunneling model, but the peak amplitude is much larger than 
expected from the available hole band parameters. Data as a function of parallel magnetic field 
reveal additional anomalous features, including a recurrence of a zero-bias tunneling peak at very 
large fields. In a perpendicular magnetic field, we observe a robust and narrow tunneling peak at 
total filling factor vt = 1, signaling the formation of a bilayer quantum Hall ferromagnet. 

PACS numbers: 73.43.Jn, 73.40.Ty, 73.21.-b 



Measurements of tunneling between a pair of 
low-disorder, GaAs two-dimensional electron systems 
(2DESs) have revealed some of the most basic [3, 12, 0, 0, 
0,0], as well as exotic [3] properties of these systems. The 
measurements have probed, e.g., the Fermi contour of the 
2D electrons ll| , and have determined their quantum life- 
time 0, i, S, E i, i, 0, [HI, El . More recently, tunneling 
experiments on closely spaced, interacting electron bilay- 
ers at the total Landau level filling factor vt = 1 revealed 
an unusual, excitonic quantum Hall ground state [7j . 

Here we report an experimental study of the differen- 
tial tunneling conductance between very dilute, closely 
spaced GaAs 2D hole systems (2DHS) as a function of 
interlayer bias [HI]. In each of the layers, the inter-hole 
distance measured in units of the effective Bohr radius, 
r s = -. — m J; ; — , reaches up to ~ 13, attesting to the 

extreme diluteness of the 2DHS (here, p is the 2D hole 
density in each layer, to* = 0.2too is the hole effective 
mass [1J], and e = 12.4 is the GaAs dielectric constant). 
The physics of such a dilute system is expected to be 
dominated by interaction effects. Indeed, in the pres- 
ence of a perpendicular magnetic field (B±) and at total 
Landau level filling factor Vt = 1, we observe a strong 
zero-bias tunneling peak which persists up to a temper- 
ature of T ~ 0.6K, indicative of significant interlayer 
and intralayer interaction in this system 0] • Despite the 
extreme diluteness of our bilayer 2DHS, the shape and 
temperature-dependence of the tunneling peak in zero 
applied magnetic field are consistent with an existing 
model based on Fermi liquid theory [1, 0]. The peak 
amplitude, however, is roughly 8 orders of magnitude 
larger than that calucated using the conventional hole 
band parameters. Our measurements in samples with 
different barrier thicknesses suggest that the hole mass 
relevant to tunneling through an AlAs barrier is much 
smaller than expected. Finally, tunneling measurements 
performed with an applied parallel magnetic field B\ | re- 
veal another puzzle: while the amplitude of the zero-bias 
tunneling peak quickly decays with B\\ consistent with 
the expected momentum conservation requirement, the 
peak reappears at very large B\ | . 



We performed interlayer tunneling measurements on 
double quantum well samples grown on GaAs (311)A 
substrates. The structures consist of a pair of 15 nm- 
wide GaAs quantum wells, seperated by 7.5, 8.0, or 11 
nm AlAs barriers, and flanked by Si-modulation-doped 
layers of Alo.21Gao.79As. While we will discuss the data 
from all three samples, the data shown in this report are 
all from the sample with a 7.5 nm barrier. The as-grown 
density for this sample is ~ 2 x 10 10 per layer, and the 
typical mobility at 0.3 K is ~ 20 m 2 /Vs. A Y-shaped 
mesa, illustrated schematically in the inset to Fig. 1(a), 
with an active region of 100 /.im x 700 /im, was etched 
into the wafer. Alloyed InZn contacts at the ends of the 
seven arms allow for making independent contacts to in- 
dividual 2D hole layers using the selective gate depletion 
scheme (TBI ] . This Y-configuration enables us to measure 
electrical transport in each of the two layers seperately, 
and to balance their densities using front and back gates 
which cover the active area, before proceeding with the 
interlayer tunneling measurements. The latter were per- 
formed by applying a low-frequency, 18/iV peak-to-peak 
ac excitation on top of an interlayer dc bias to the top 
layer, and measuring the current in the back layer using 
standard lock-in techniques. The lock-in phase was set 
by confirming that the measured differential (ac) conduc- 
tance matched the derivative of the dc I — V curve for 
the tunneling measurement. 

As illustrated in Fig. 1(a), the differential conduc- 
tance of the tunnel junction formed by the pair of 2D 
hole sheets with equal carrier density shows a pronounced 
resonance near zero interlayer bias. This zero-bias reso- 
nance is remarkably similar to that seen previously in 
bilayer electron samples @, 0, 0, 0, 0] , despite here being 
taken in the extremely dilute limit (r s ~ 11 — 13 for holes, 
r s ~ 1 — 6 for electrons). In the earlier works [3,0,0], the 
dc conductance, I/V, derived from the data was found 
to closely match a theoretically derived [1, 0] expression: 
{?{V) — St 2 ^^? ( eV y i+r 2 , where S is the area of the 
tunnel junction, t is the tunneling matrix element, to* is 
the density-of-states effective mass for the 2D carriers, T 
is inversely proportional to the quasiparticle scattering 
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FIG. 1: (a) The differential conductance of the bilayer hole system as a function of interlayer bias. The layers have equal 
densities of p = 2.4 x 10 10 cm -2 , corresponding to an r s of 11, and are seperated by a 7.5 nm thick AlAs barrier. The solid 
line represents a fit of Eq. (1) to the data. The inset shows a schematic diagram of our Hall bar. Front and back gates cover 
the active (gray) area to allow for tuning the density of each layer. Contacts (crosses) can be made to either the front (F) or 
back (B) layer by energizing selective, depletion (contact) gates (not shown) on the back or front of the sample, respectively. 
Tunneling measurements were performed by energizing all the contact gates, and using only the four contacts in the middle of 
the Hall bar. (b & c) Tunneling amplitude (A) and width (F), derived from fitting the data to Eq. (1), are shown as a function 
of temperature for three layer densities: p = 1.8 X 10 10 cm~ 2 (r a = 13), 2.1 x 10 10 cm" 2 (r s = 12), and 2.4 x 10 10 cm" 2 (r s = 11). 
Definition of symbols is provided in the inset to (d). (d) The width T, scaled to the Fermi energy (ef), is shown as a function 
of the temperature T, scaled to the Fermi temperature (Tp). The solid line is the theoretical prediction from Fig. 8 of Ref. 0, 
offset vertically to account for the disorder scattering in our sample at p — 2.4 x 10 10 cm -2 . 



lifetime, and both T and \eV\ are assumed to be much 
smaller than the Fermi energy Ep. This zero-bias reso- 
nance results from a close matching of the momenta and 
energies of the carriers in the two wells. In Fig. 1(a), 
we show a fit of our data to the differential conductance 
derived from this expression (ljj . 



dV [ ' [r 2 + (eT/) 2 ] 2 



(1) 



where A = St z 'f . 

n Tin 



2 2em_ geen j n pjg i( a ) ; this expres- 
sion provides a remarkably good fit to the experimental 
data, despite a T (0.15 mV at 0.31 K) and applied bias 
which are comparable to the Fermi energy (0.29 mV). 

In Fig. l(b & c) we summarize the T-dependence of 
values of A and T extracted from the fits of our data to 
Eq. (1). In general, the width r is inversely proportional 
to the quasiparticle scattering lifetime. As seen in Fig. 
1 (c) , the fitted widths for all three carrier concentrations 
appear to extrapolate to finite values as T — > 0. Static 
disorder, mainly from ionized impurities, is believed to 
lead to a finite zero-T scattering lifetime, and hence a 
finite tunneling peak width 0, [1| . Our observed decrease 
in the extrapolated zero-T width with increasing carrier 
concentration is likely a result of the improved screening 
of the disorder, and is consistent with previous bilayer 
electron data. The increase of T with increasing T in 
the bilayer electron case has been attributed to carrier- 
carrier scattering [1, 0, 01- Jungwirth and MacDonald 
j9j calculated this scattering, and found excellent, quan- 



titative agreement with the electron bilayer data [4j . As 
shown in Fig. 1(d), this prediction underestimates the 
change in the measured width as a function of tempera- 
ture for our data. The discrepancy is possibly due to the 
very large r s in our system [9j. 

In Fig. 1(b) we observe that the amplitude A is T 
independent to within ±2.5% for the three densities in- 
vestigated here. According to Eq. (1), the amplitude 
A is proportional to the square of the tunneling matrix 
element t, which is intuitively T independent. In a single- 
particle picture, t is the hopping matrix element between 
the two planes, and is related to the energy splitting be- 
tween symmetric and antisymmetric wavefunctions span- 
ning the two wells, 2t = As as- Using a density-of-states 
effective mass m* = 0.2mo 1J| , we find Ag^s = 8.6 nV 



for p = 2.4 x 10 cm and slightly smaller for lower 
p (see Fig. 1(b)). These Ag^g values are 4 orders of 
magnitude larger than expected from simple estimates of 
Ag^s based on the structure of our bilayer system and 
the hole band parameters 171 ]. 

Are strong interlayer correlations responsible for the 
anomalously large zero-bias tunneling peak that we ob- 
serve in our very dilute bilayer sample? Note that in our 
sample the ratio of the average in-plane interparticle dis- 
tance (r) = l/*Jp to the interlayer distance d is about 
2.9-3.3 (we define d as the center-to-center distance be- 
tween the two quantum wells). If the large tunneling 
peak is indeed related to the diluteness of our bilayer 
layer system, then we would expect a similarly large peak 
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FIG. 2: (a) (Color online) False color plot of the differential 
tunnel conductance of the hole bilayer as a function of in- 
terlayer bias and applied parallel magnetic field. Data were 
taken at 0.31K and a hole density of 2.1 x 10 10 cm -2 per layer. 
For clarity, the tunneling traces shown in (b) have been offset 
vertically (in steps of 2/iS). 



in samples that have a similar (r) jd ratio but whose bar- 
rier is slightly thicker. Our measurements on two other 
samples with (r)/d ~ 3, however, indicate otherwise: In 
a sample with a barrier width u>b = 8.0 nm we observe 
a tunneling peak whose amplitude is a factor of ~ 10 
smaller than the peak for the 7.5 nm-barrier sample, and 
in a sample with % — 11 nm there is no measurable 
tunneling. These observations are, on the other hand, 
consistent with a single-particle picture in which A 949 
is exponentially sensitive to the width of the barrier [1 81 ] . 
To be quantitative, using a self-consistent Shroedinger 
solver for our sample structure, we find that the calcu- 
lated values of As^s approach those extracted from the 
measured data for both 7.5 nm- and 8.0 nm-barrier sam- 
ples if we take m z — 0.19too as the mass in the growth 
direction, instead of the expected value m z = 0.63mo 
While our assumption of a single mass for the 
growth direction means we cannot directly interpret this 
value as being the mass of holes in AlAs in the growth di- 
rection, our data do strongly suggest that the hole mass 
relevant to tunneling through an AlAs barrier is smaller 
than has been suggested in the literature [17] . 

The zero-bias tunneling peak arises from a close match- 
ing of both the energy and momentum of the carriers in 
the layers. By applying a magnetic field B\\ parallel to 
the layers, the in-plane (canonical) momenta of the car- 
riers in one well shift with respect to the other by k& = 
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FIG. 3: (Color online) The differential tunnel conductance 
as a function of interlayer bias is shown at a layer density 
of 2.1 x 10 10 cm" 2 in a perpendicular field of S x =1.72 T, 
corresponding to a total (bilayer) filling factor of vt = 1. 



(eB\\d)/h [20l| . We would expect the zero-bias peak to 
have maximum amplitude at zero field, where the Fermi 
contours of the two layers overlap, and decrease rapidly as 
the carrier momenta in one layer are shifted with Bit and 
no longer overlap the momenta in the other layer l|, @] . 
As shown in Fig. 2, we find that there is indeed a rapid 
suppression of the zero-bias peak upon the application of 
a small B\\ [2l[, of the order of 1 T. However, Fig. 2 re- 
veals a very rich and puzzling evolution of the tunneling 
spectra at higher B\\. If we assume parabolic in- plane dis- 
persions, the zero-bias tunneling peak should disappear 
beyond B?, = (2hk F )/(ed) = 2. IT, where k F = v / 27rp 
is the Fermi wavevector assuming a spin-degenerate sys- 
tem [22[. However, well beyond this field, we find peaks 
which evolve smoothly as a function of bias and i?| | , with 
two dispersing peaks intersecting near zero bias at ~ 9T 
and ~ 14T. This suggests a strong distortion of the hole 
bands and Fermi contours in a large B\\. In particular, 
a simple change in the band structure curvature, due to 
the finite hole layer thickness [23j , would not explain the 
re-emergence of a zero-bias peak at high Bu. 

Finally we present tunneling data in our bilayer hole 
system in a perpendicular magnetic field, B±. As has 
been shown in bilayer electron samples with strong inter- 
layer interaction, when B± is tuned such that the total 
Landau level filling factor {vt) is one, a many-body zero- 
bias tunneling resonance ensues Q ■ With each layer at a 
filling factor of 1/2, if the ratio of the in-plane to the out- 
of-planc Coulomb energy, parameterized by d/ls (where 
Ib is the magnetic length) is < 1.8, the occupied states 
in one layer are spatially correlated with vacant states 
in the other layer, and tunneling is resonantly enhanced 
0| • This pairing has been further confirmed in both elec- 
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tron and hole bilayers using counterflow transport stud- 



ies 124, 25 



261 127| . Here we provide tunneling data (Fig. 
3) for this unique state for our bilayer 2DHS. Despite 
a higher background level than the zero-field tunneling 
data, there is a strong zero-bias tunneling peak, persist- 
ing up to T ~ 0.6 K. This is higher than the temperatures 
where the equivalent peak in bilayer electron samples dis- 
appears 0, [28| and reflects the very strong pairing in our 
sample: the parameter d/ls = 1.1 for our sample is in- 
deed smaller than cI/Ib > 1-5 for the published electron 
bilayer data. This persistence of the tunneling peak up to 
higher temperatures is also consistent with the counter- 
flow measurements in bilayer hole samples with relatively 
small d/ls values [H,!!?]]. We defer a detailed discussion 
of the vt = 1 bilayer hole tunneling peak height and 
width on temperature and d/ls to a future communica- 
tion. 

In conclusion, we have demonstrated that certain as- 
pects of interlayer tunneling in dilute bilayer 2D hole 
systems are remarkably similar to what has been seen 
in bilayer electrons. The shape and width of the zero- 
field tunneling peak can be well explained in terms of 
a Fermi liquid theory developed to understand the elec- 
tron data. The anomalously large size of the zero-field 
peak is consistent with a smaller than expected hole mass 
relevant to tunneling through an AlAs barrier. The par- 
allel magnetic field dependence of the tunneling spectra 
reveals surprises when compared to what is qualitatively 
expected based on a simple band structure. Finally, data 
taken in perpendicular field show enhanced tunneling at 
i>t = 1 in bilayer holes at relatively high temperatures, 
consistent with the system being deep in the strongly- 
interacting regime. 

We thank the Princeton NSF MRSEC and DOE for 
support, and D.A. Huse and R. Winkler for illuminating 
discussions. 
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